Coexisting orders in the quarter-filled Hubbard chain with elastic deformations 
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The electronic properties of the quarter-filled extended Peierls-Holstein-Hubbard model that 
includes lattice distortions and molecular deformations are investigated theoretically using the 
bosonization approach. We predict the existence of a wide variety of charge-elastic phases depending 
of the values of the Peierls and Holstein couplings. We include the effect of the Peierls deformation 
in the nearest-neighbor repulsion V, that may be present in real materials where Coulomb interac- 
tions depend strongly on the distance, and we show that the phase diagram changes substantially 
for large V when this term is taken into account. 
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I. INTRODUCTION 

The electronic properties of the low dimensional in- 
teracting electron systems have attracted great interest 
for many different reasons, the main being that they are 
simpler to analyze than higher dimensional ones and they 
could then be used to gain some insight on their higher 
dimensional counterparts. Typically this connection can 
be envisaged by coupling many ID systems to build up 
a higher dimensional array. They are also interesting 
on their own since many materials show a quasi-one- 
dimensional behavior in a certain range of temperatures 
where higher dimensional couplings can be neglected. A 
third motivation to study ID systems comes from the 
very exciting cold atom systems which provide a test lab- 
oratory for many important theoretical developments al- 
ready done and also motivates further analysis that could 
be in principle directly tested in real systems. 
In particular, and in connection with real materials, 
quarter-filled systems have received a lot of attention 
since they are very good candidates to describe the so 
called organics charge transfer salts like the Bechgaard 
salts (TMTSF) 2 X or (TMTTF) 2 A whre X=PF 6 , AsFgk 
The Hubbard model is the basic model for understand- 
ing the electronic properties of these quasi 1-D organic 
conductors and in some cases, lattice deformations could 
play a central role. In addition to this, in crystalline 
materials where the building block of the crystal struc- 
ture is a large molecule, the vibrational properties of the 
molecules often have large effects on the electronic prop- 
erties of the material^. These materials exhibit a variety 
of electronic states with spatially inhomogeneous charge, 
spin and lattice structures. 

Several experimental studies have shown clear evidence 
for the existence of charge oidev^^^^^- in organic 
charge transfer solids and this has stimulated consider- 
able theoretical effortsii^ 2 -^^^^^^. These mate- 
rials involve Coulomb repulsion, both on-site and be- 
tween nearest-neighbors (U, V respectively) and arc 1/4- 



filled. Due to the nearest-neighbors Coulomb repulsion, 
it is expected that in these materials a charge ordering 
• — o — • — o sets in, where • and o correspond to occu- 
pied and unoccupied sites. The charge order • — o — • — o 
corresponds to the so called 4kp Charge-Density- Wave 
(CDW), where Uf = Tr/4a with a the lattice constant. 
However, for large values of V, and when electron-phonon 
interactions arc included, the ground state corresponds 



to an ordering like • — • — o - 



,11.12.13.14.15.16 



The charge 



order corresponds to a tetramerized 2kp CDW and it co- 
exists with the 2kp periodic modulation of the intersite 
distances, called Bond-Ordered- Wave (BOW) or a mixed 
2kp + Akp BOW. This charge order is called Bond- 
Charge Density Wave (BCDW) . On the other hand, if the 
Coulomb interaction between two sites r, and rj , is expo- 
nentially short ranged, i.e. of the form exp[—o:|ri— rj |]/r^ 
with cv,£ > 0, lattice distortions may have an impor- 
tant effect on the Coulomb interactions between nearest- 
neighbor. 

Motivated by the previous discussion we study, using 
the Abelian bosonization approach and a semiclassical 
analysis of the effective theory, the ground state phase 
diagram of the Peierls-Holstein-Hubbard model includ- 
ing the effects of the Peierls deformation in the nearest- 
neighbor Coulomb interaction V. We show that the 
ground state at zero temperature indeed favors several 
period two and period four distortion patterns, stem- 
ming from a competition between elastic and charge- 
spin energy. These patterns spontaneously break transla- 
tional symmetry, with different phases depending on the 
Coulomb interaction strengths U and V and on the value 
of the spin-phonon couplings. A rich phase diagram is ob- 
tained, including all the combinations of Charge Order 
and 2k p — Akp deformations function of the spin- 
phonon couplings. 

When the effect of the Peierls deformation in the 
nearest-neighbor Coulomb interaction is not considered, 
our analysis reproduces previous results obtained with 
numerical methods^ . The main novelty in our treat- 
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ment is the inclusion of a nearest neighbors Coulomb 
repulsion V which depends on the distance between ad- 
jacent ions and hence depends on the lattice distortions. 
The phase diagram then changes substantially indicating 
that this variation on V should not be neglected in com- 
paring to experiments. Our analysis allows for a descrip- 
tion only at a qualitative level, but it has the advantage of 
being easily gcneralizable to include other perturbations 
as longer range Coulomb interactions, to treat arbitrary 
filling fractions, to add an external magnetic field, etc. 

The paper is organized as follows. In Section II we 
present the model and its analytical treatment. The 
phonon sector is described in the adiabatic approxima- 
tion by classical static deformations and we analyze the 
possible patterns. In Section III, we use the bosonization 
approach to represent the electrons operators in terms of 
bosons and analyze the effective description by consid- 
ering all the relevant perturbation terms as semiclassical 
potentials. We draw a qualitative phase diagram with 
from analisys. Special emphasis is put on the character- 
ization of the ground state phases that result from the 
combination of charge-spin degrees of freedom and clastic 
effects in different parameter ranges. Finally, in Section 
IV we present a summary and conclusions of the present 
work. 



II. MODEL HAMILTONIAN 

We consider a ID extended Hubbard model with 
Hamiltonian given by 

3,0 

+ U X! n 3>l n 3,l + V 3^ + 1 n 3 n 3+l + 

3 3 

3 

where Cj j(T denotes the annihilation operator of an elec- 
tron with spin cr(=t, I) at the jth site, and rij = n^f + 
rijj is the charge density with rij !<7 = c] a Cj a . £i,i+i is the 
hopping integral, the parameters U(> 0) and Vj 7+1 (> 0) 
denote the magnitudes of the on-site and nearest neigh- 
bor interactions and j3j represents the internal molecular 
deformation. We focus on the ^-filled case, with the av- 
erage number of electrons per site (rij) = 1/2. 

The interaction of the charge degrees of freedom in a 
homogeneous Hubbard chain {tjj+i = t) with phonons 
is usually modeled by a linear expansion of the exchange 
couplings around the non distorted values t, 

tj.j+i — *(1 _ 9p Spj) 

where 5pj is the local lattice distortion and gp is the 
electron-phonon interaction. On the other hand, local 
deformations of the molecules can produce changes of the 



on-site (or molecular) orbital energies and can simply be 
taken into account by a Holstein term 

9H^2S H ,jrij, (2) 
3 

where gn is the electron-phonon coupling constant of the 
on-site type (in Eq. (fT]) we change f3j — ► gn $H,j)- 

In addition to the purely electronic terms in Eq. (TTJ) 
and unlike previous analysis, we include the effects of 
the Pcicrls deformation in the nearest-neighbor Coulomb 
interaction V that may be important when the Coulomb 
interaction between two sites r, and Vj is exponentially 
short ranged, i.e. exp[— a\ri — Tj\]/r^. Specifically, we 
perform a linear expansion of the Vjj+i couplings around 
the non distorted values V, 

v 3,j+i - V(l-gvS P ,j) 

where gy is the coupling constant which measures the ef- 
fects of the bond deformation on the Coulomb repulsion. 
The complete Hamiltonian, including the elastic energy 
in the adiabatic approximation, is written as 

H = - f J2( l - 9P 8p,3)[Cj >tr Cj+l,<T + c)+i,a c jA + 
3,v 

+ U n 3A n 3-i + V ~ 9V n 3 n 3 + ^ + 

3 3 

■ !IU X Sh -3 U 3 
3 

+fE^ + TE^ (3) 

3 3 

where the parameters Kh and Kp are the conventional 
elastic constants for the on site deformation and the lat- 
tice distortion, respectively. 

Now, we consider the lattice deformations. At quarter- 
filled ((rii) = 1/2) the most general period four lat- 
tice deformations without collective displacement, can be 
parametrized as 

S pj = Sp y d cos(irxj) + Sp t t cos(7r/2 Xj + £ — 7r/4) 
Shj = S H ,d cos(tt xj) + S H ,tCos(Tr/2xj - C) (4) 

where Xj is the position of the j — th site (lattice con- 
stant a has been set to one), and the phases £ and Q de- 
termine the spatial patterns of tetramerization Spj and 
5ji,t respectively, while 6p t d and 5h,<i are the amplitude 
of the dimerization in the lattice distortion and the intra- 
site electron-phonon deformation. This parametrization 
is the most general supported by bosonization, as such 
a deformation is commensurate. Period four deforma- 
tions cause commensurability of relevant perturbations 
at (rij) = 1/2 and provide a mechanism for a charge gap 
in this regime^. Numerical evidence of the dominance of 
period four lattice deformations has been obtained from 
self consistent computations^. A uniform deformation, 
leading to global size change, can also appear; this would 
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produce a uniform shift in t and V, which is inessential 
to our present analysis. 

Before analizing the case gv > 0, we summarize all previ- 
ous numerical results in the 1/4 filled case with gv = 0. 
When gp = gn = and U — ► oo there exists a crit- 
ical value of V (V c = 2t) for the appearance of a 4k f 
Charge-Density- Wave if V > V c Fig. Q^c). For V < V c 
the extended Hubbard model corresponds to a Luttinger 
liquid (LL) with no charge order. For finite U, the value 
of V c has been calculated withing strong coupling pertur- 
bation theory 2 ^ and numerical methods 2 ^. 
For small positive gp and gn, and V < Vc, the ground 
state has a dimerized AkF Bond-Order- Wave (BOW) 
with uniform site charges Fig. QJa). If gp is sufficiently 
large (gu fixed), there is a superposition of the 2kp and 
Akp BOWs accompanied by the • — • — o — o charge or- 
der. This state correspond to the Bond-Chargc-Dcnsity- 
Wave state (BCDW) Fig. []Jb). On the other hand, at 
large gn and small gp, a Akp CDW is found, with uni- 
form bonds Fig.QJc). If we increase gp, the ground state 
becomes the Akp CDW-2fcp SP, where the charge order 
• — • — o — o, is accompanied by a 2kp lattice distortion 
Fig. 0Jd). For V > Vc, the Akp BOW phase is dimin- 
ished due to the strong Coulomb repulsion V and phases 
with charge-order dominate the phase diagram. In this 
work we present new results for a 1/4 filled band in- 
cluding the electron-phonon interaction (Peierls and Hol- 
stein couplings) and lattice deformations in the nearest- 
neighbor Coulomb interaction V that has not previously 
been studied. We show that this term may change sub- 
stantially the phase diagram for large V. 
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FIG. 1: Possible lattice deformations and charge density wave 
ground state in the Peierls-Holstein-Hubbard model at quar- 
ter filling. The Grey, black and white circles correspond to 
site charges of (n) = 0.5, 1 and 0, respectively. We reproduce 
the previous results of^ and we show two new phases (e), (/) 
that appear when gv > 0. In this figure (a) The AkF BOW 
state, with dimerized bond orders and uniform site charge 
densities, in (b) represent the BCDW state, with bond orders 
2Uf > AkF and charge order • — • — o — o. (c) The 4/cf CDW 
with • — o — • — o charge order and uniform bond order, (d) 
The AkF CDW-2fcp SP state, with • — o — • — o charge order 
and bond orders 2&f. For gv = 0.2, we have phase (e) i.e., 
AkF CDW-4fcp SP with • — o — • — o charge order together 
with dimerization bond order. Finally, in the phase (f) (AkF 
CDW-4fc F SP + 2k F CDW-2fc F SP) is characterized by three 
different site charges | — e' , i + e, | — e', \ + e" with bond 
orders 2/cf > AkF- 



III. BOSONIZATION APPROACH AND 
SEMICLASSICAL ANALYSIS 

In order to analize semi-quantitatively the low energy 
properties of the model given by Eq. ([3]), we use the 
Abelian bosonization method which is generally powerful 
for the description of one-dimensional chains (for further 
details see Ref i 20 ' 21 i 24 ' 25 i 26 ). Dimensionless parameters, 
are used below. They are introduced using t as the en- 
ergy scale as follows: gp — > gp = (t/ Kp) x / 2 qp, gn — > 
9h = {t/K H ) 1/2 g H , gv -> gv = {t/Kpf^gv, Sp, s 
~5 P , S = (Kp/ty/ 2 S P , a and 5 H , S -> ~5 H , S = [K H /t) l ' 2 5 H ,s 
(,s = d,t). To obtain the corresponding low-energy, we 
write the fermion operator as 2 ^ 



_ e ikF,„X e -iV4TT<f>L,cr(x) 

_\_ e -ikF,crX e iV4lr<f>R }tr (x) _j_ ^ 



where kp /j are the Fermi momenta for up and down spin 
electrons and (j)R,L,iT are the chiral components of two 
bosonic fields, introduced as usual in order to bosonize 
the spin up and down chiral fermion operators ipRL,a- 
The dots stand for higher order terms 2 ^. They take 
into account the corrections arising from the curvature of 



the dispersion relation due to the Coulomb interaction. 
For non-zero Hubbard repulsion U and V , the low en- 
ergy effective Hamiltonian corresponding to (|3|) written 
in terms of the bosonic fields 0| and 4*i has a complicated 
form, mixing up and down degrees of freedom. We define 
4>cr = 4>R,a + ^L.cr and introduce linear conbinations of 
to describe the charge and spin degrees of freedom, 
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(<f>l +cj>i) , (j> s 



1 



Or _ <t>i) , 



Then, we can rewrite H = dxJi as (at h = 0): 

H = Hph + Hf re e + V e f f 



(6) 



(7) 



where 

T~tph 

Hfree 



Veff = 



1(1 

Vc_ 

An 
Air 

V c + 



2 y 2 



K, 



\d x (j) c 



K p {d x e c ) 



1 



K, 

V s + Vi n t 



[d x( t> s ) 2 + K s {d x e s f 



(8) 
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V c = ffi/4 cos 40c + Ai gp 6 P ,d sin 2<j> c + A 2 <?# 5 Htd cos 20 c + A 3 g v <5p,d cos 20 c 
V ff = g s cos20 s 

Vi„ t = —gp 5p,t cos (0 C - £) cos S - <5f/, t sin (0 C + £) cos (f> s - A 4 5 P ^ cos 20 c cos 20 s 

+A 5 g v 8p tt cos S sin ( 3 4> c + £ - — J - A 6 gy £p,t cos S sin ( <p c — £ + — J (9) 



plus several vertex operators that are kept only when 
they are commensurate (non-oscillating in space) and 
constitute relevant perturbations to the Gaussian con- 
formal field theory. In Eqs. §E§ and © <f> c and </> s are 
the charge and spin fields respectively, (9 CjS are the dual 
fields defined by d x 9 CiS = dt4>c,s and the phases £ and 
C determine the spatial patterns of the tetramerization 
(see Eq. (|3|)). The parameters v c and v s are the ve- 
locity of the charge and spin excitations, and K c and 
K s are the corresponding Tomonaga-Luttinger parame- 
ters. From the perturbative calculation presented in^i 
the coupling constants are given by giu oc U 2 (U — 4 V) 
and g s oc U — aU(U — 2 V) where a is some numerical 
constant. Besides, 



g yA oc U 2 (U~4V) 

g s oc U - a U(U - 2 V) 

Ai a U 

A 2 oc U 

A 3 oc (U-cte)V 

A4 oc V 

A 5 oc UV 

A 6 oc (U-cte)V. 



per site depends on four configuration parameters, 
e(0 C) 8 ,£,C) = Veff{<f>c,<f>s,ti,0 



(11) 



The Eq. ([9]) is the bosonic self-interaction poten- 
tial defining a Multi-sine Gordon theory. Extensive 
analysis of the competition between different harmon- 
ics in multi-frequency sine-Gordon theories have been 
performed i n 28 i 29 i 30 , mainly focused on the double sine- 
Gordon model. The three-frequency case has also been 
recently discussed i n 31 i 32 i 33 . For our purposes it will be 
enough to perform a scmiclassical treatment, as detailed 
in the next Section. 

In order to perform the semiclassical analysis of the 
ground state and present a schematic phase diagram, we 
use the bare coefficients in Eq. ([9]) and we assume the 
qualitative phenomenological dependence of the coupling 
constants <?i/4, g s and A's on the microscopic parameters 
as infTUl 

The aim of the present work is to search for the possi- 
bility of elastic deformations that lower the energy with 
respect to the homogeneous non-deformed case. The sim- 
plest analysis of the effective theory Eq. 0, ([5]) and 
which has proven to be useful in related case o 32 ' 33 ! 34 ' 35 , 
consists in treating the self-interaction terms in Eq. @ 
as a classical potential to be evaluated in constant field 
configurations. Within this approximation the energy 



so that the minima can be found analytically. In the pre- 
vious expression, c6 c and <j) s are the charge and spin fields, 
and the phases £ and £ determine the spatial patterns of 
the tetramerization. 

Our results on coupled Holstein-Peierls-Hubbard chain 
are summarized in Fig. [3J First, we consider the case 
gv = 0. From the perturbative calculation, the coupling 
giu is positive for small V but becomes negative for large 
V. When g X j A > 0, the phase variable C is to be fixed 
at <p c = 7r/4, while for gin < 0, <f> c = 0. 

Let us consider finite values of gp and gn- In Fig. 
[3{a) we have set 51/4 > (or U > 4 V), and we increase 
gp from zero (for fixed small cjh)- Unlike what has been 
found numerically ini2 (where no distortion patterns were 
favored, probably due to finite size effects), the phases 
(4> c ,4>s) are locked in (f-,f) (independently of £,C)- In 
this case, the ground state becomes the 4k p BOW while 
the site charge density remains uniform ((rij) = 1/2), 
the bond-order is now inhomogeneous and has the form 
shown in Fig. QJa). For gp larger than a critical values 

(sf^Tr) and the 



(10) 9p ^,9p(§h), wehave(£,Cj 



BCDW depicted in Fig. [TJb) becomes the ground state. 
Hence the bond-order pattern now is t(l + a), til — a'), 
t(l + a), t(l + a") and the charge order follows the o — 

• — • — o pattern. Now consider the situation when gn 
reaches a critical value gji(gp)- In that case, (cj) c ,4> s ) = 
(0, f) (independently of and the Ak F CDW with 

• — o— •— o charge order and uniform bond order becomes 
the ground state Fig. [TJc). The SP distortion, which 
corresponds to the 4k F GDW-2kp SP state (shown in 
Fig-HUd) occurs only for large values of gn and gp w 0.8. 
In this case (€X,<l>p,M = (f , 7r, 0, 7r). 

In Fig. [3fb) we can see several important trends when 
ffi/4 < (or V > C//4). First, the size of the 4fc F BOW 
region decreases considerably with decreasing g-i/ A < 0, 
showing that it takes stronger e-ph coupling to overcome 
the "natural" tendency towards • — o — • — o charge 
order for small gp. For the same reason, the sizes of 
the Ak F CDW and Ak F CDW-2/cf SP regions increase 
considerably when g\/ A is negative. Next we examine 
the effect of gv When gv = 0.2 and g^u < 0, we 
can see that the 4kp BOW region "disappears" and for 
9p < I- 5 the new phases Ak F CDW-4fc f SP and Ak F 
CTm-Ak F SP + 2k F CDW-2fc F SP dominate. Unlike 
the previous case with gv = 0, in the these phases, there 
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FIG. 2: Semiclassical energies for the minima in terms of gp 
for gi/i — —1.5, gu = 2.2 and gv = 0. 
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IV. SUMMARY AND DISCUSSIONS 

We have investigated the electronic ground-state 
phase diagram of the one-dimensional Hubbard model at 
quarter filling with on-site (U) and nearest-neighbor (V) 
Coulomb repulsion coupled to adiabatic phonons that 
includes Holstein and Peierls interactions. In addition to 
the conventional terms, we have included the effects of 
the Peierls deformation in the nearest-neighbor Coulomb 
interaction V that may be important when the Coulomb 
interaction between two sites and rj is exponentially 
short ranged, i.e. cxp[— a\i*i — rj\]/r^ > 0). We 

have used Abclian bosonization to analyze the low 
energy properties of the model and we have performed 
a semiclassical analysis of the bosonized effective theory 
to identify the possible phases of the ground state. We 
have found that several charge orders combined with 
lattice deformations describe the ground states of the 
system depending on the microscopic parameters U, V 
and the spin-phonon couplings gp, c/h and gv In each of 
these phases a non trivial elastic deformation is favored, 
grouping together blocks of two or four sites, while the 
charge sector adopts different charge orders. 
A detailed analysis for particular values of gv, shows 
the following phases at zero temperature: 
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FIG. 3: Schematic phase diagram in the plane gp vs §h for 
different values of gv- (a) gi/4, = 1, gv = 0; (6) gi/ 4 = —3, 
gv = 0; (d) gi/4 = —3, gv ~ 0.2. The boundaries of the 
phases result from level crossings as shown in the Fig. 
The case with §1/4 > and gv > is similar to the case (a) 
and we don't show here. The effect of the Peierls deformation 
in the nearest-neighbor Coulomb repulsion is more important 
for large V, that is, §1/4 < . 



is now a small dimcrization in the bonds (Fig. [He) and 

Fig.rnf)). 



< For U > V (31/4 > 0), gv = and for small 
and nonzero gp and gn, the ground state has 
a dimerized Akp Bond-Order-wave (BOW) with 
uniform site charges Fig. [Tfa). If is sufficiently 
large (gn fixed), there is a second dimcrization 
of the dimerized lattice, leading to the BCDW 
state (a superposition of the 2fcp and 4/cf BOWs), 
accompanied by the • — • — o — o charge order Fig. 
QJb). On the other hand, at large gn and small 
gp, a Akp CDW is found, with uniform bonds Fig. 
[ljc). If we increase gp, the ground state becomes 
the Akp CDW-2fcf SP, where the charge order 
• — • — o — o, is accompanied by a 2kp lattice 
distortion Fig. [TJd). 



For V > U (g 1/4 < 0), the 4fc F BOW phase is 
diminished due to the strong Coulomb repulsion V 
and phases with charge-order dominate the phase 
diagram Fig. ©. 



• For V > U (<7i /4 < 0), g v > 0, the 4k F BOW dis- 
appears because the Coulomb interaction between 
nearest-neighbor is increased due the Peierls defor- 
mation. The CDW's phases has a small contribu- 
tion of Akp BOWs combined with the original 2kp, 
giving rise to two new phases: Ak F CDW-4fcp SP 
and 4k F CDW-4fcp SP + 2k F CDW-2fc f SP 

In conclusion, we have argued that the coexistence of 
charge order and lattice distortions is dominated by Akp 
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BOW, BCDW for g 1/4 > 0. On the one hand, when 
5i /4 < 0, the 4fc F CDW and 4fc F CDW-2fc F SP dominate 
if 9p < 9 c p, Fig- Hb); for g P > g c p the BCDW phase 
dominate. On the other hand, when the phonon-latticc 
coupling, gv is taken into account (only for g\u < 
because is necessary V > U), two new phases that we 
named 4fc F CDW-4fc F SP and 4fc F CDW-4fc F SP+2fc F 
CDW-2fc F SP (Fig. [He) and Fig. [2(f) respectively) dom- 
inate the phase diagram. These phases differ from the 
original 4fc F CDW and 4fc F CDW-SP, in a small dimer- 
ization in the bonds, Sp,d and a small contribution of 2kp 
CDW in the last one. 

Numerical analysis of the present system including 
Peierls deformations in the nearest-neighbor coulomb in- 
teraction would be helpful in order to determine the 
phase boundaries quantitatively. We hope that our pre- 
dictions will stimulate further X-ray or neutron scat- 
tering experiments on compounds like Bechgaard salts. 
They should be able to distinguish the different distortion 



patters predicted here. Since the effects of the Peirels de- 
formations in the nearest-neighbor coulomb interaction is 
proportional to V (see Eq. ((T|)), we expect that the ef- 
fects of this term will be more important in compounds 
where V is larger. 
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